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ABSTRACT 

We performed a detailed analysis of the use of [C/N] measured in red giant branch stars between the completion of the first dredge 
up and the red giant branch bump ([C/NjpDu) as age indicator. [C/N]fdu cannot give accurate ages for individual stars, but may 
provide a general chronology for the formation of composite populations and add constraints to analyses of red giants from surface 
gravity-effective temperature diagrams. We provide a theoretical calibration of [C/Njpuu in terms of total metallicity [M/H] and age, 
for ages greater than ~ 1 Gyr, which we tested against variations in the initial heavy element distribution (scaled-solar vs a-enhanced), 
efficiency of overshooting from MS convective cores and from the convective envelopes, variations in the initial He abundance and in 
the mixing length parameter. Our calibration is compared with a small sample of available measurements of [C/N]fdu in star clusters 
and halo field stars, which at least qualitatively confirm the overall trend of the predicted [C/N]fdu with age and [M/H]. The use 
of [C/N]FDu-[M/Ft]-age relations obtained from independent sets of stellar evolution calculations cause age differences (for a given 
[C/N]fdu and [M/H] pair) up to ~2 Gyr. More accurate spectroscopic measurements of [C/N]fdu in star clusters with well-established 
ages and metallicities are required to better test theoretical calibrations of this age indicator. 
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1. Introduction 

Large spectroscopic surveys, such as LAMOST (Large sky 
Area Multi-Object fiber Spectroscopic Telescope, IZhao et ^ 
l2Q12h . SEGUE (Sloan Extension for Galactic Understanding 
and Exploration. I Yannv et alJl200^ . GAIA-ESO (iGilmore et all 
I 2 OI 2 I) . and APOGE E (Apache Poin t Observatory Galactic 
Evolution Experiment. iHoltzman et al]l20L5ll have been devised 
to provide surface gravity, effective temperature, radial veloci¬ 
ties, and photospheric abundances of several chemical elements 
for large samples of Milky Way stars, with the aim of unrav¬ 
elling our galax y detailed evolutiona r y history. In paralle l, the 
Kepler mission (iBedding et alJl201^ fchaplin et aljr2010l) pro¬ 
vides asteroseismic data for dwarfs and red giants, which lead to 
an estimation of stellar masses and radii. 

Translating this data into an evolutionary history of the 
Galaxy requires age estimates of the observed field stars. 
These are usually performed by fitting theoretical stellar evo¬ 
lution tracks to the observations in the log(g)-log(re//^) dia¬ 
gram, with additional constraints given by the measured chemi¬ 
cal abun dances and asteroseismic measurements when available 
(see, e .g., ISerenelli et al.l2013l : lRodrigues et al.l20l4lValle et all 
I 2 OI 5 I) . Ages of red giant branch (RGB) stars are particularly 
challenging because small uncertainties in ^e// (at fixed g) trans¬ 
late into large uncertainties on the star’s mass, hence its age, es¬ 
pecially when asteroseismic constraints on the mass of the object 
are not available. Additional independent constraints on the age 
of RGB stars are therefore welcome, to be used in conjunction 
with the log/gj-log/Te//) diagram (and asteroseismic masses) to 
improve the accuracy of the age estimates. 


stars: evolution - stars: interiors 


iMasseron & Gilmorel d20L5ll have recently constrained the 
age ranges spanned by thin and thick disk RGB stars using the 
[C/N] raticQ measured in RGB stars (from the APOGEE survey) 
after the completion of the first dredge-up (EDU), but fainter 
than the RGB bump luminosity, where observations (see, e.g., 
iGratton et al.ll2000l lAngelou et al.ll2015l) reveal the onset of an 
extra mixing process which continues until the tip of the RGB 
(and possibly beyond). Theoretical stellar models predict that 
this ratio -hereafter [C/N]fdu- depends on the mass of the star 
(hence its age) and to some degree on its initial chemical com¬ 
position. We present here a theoretical analysis of its calibration 
in terms of total metallicity [M/H] and age t, for t more than 
~1 Gyr, e.g. for stellar populations with well-populated RGBs, 
which host stars with an electron degenerate He-core (denoted 
here as low-mass stars). 

The next section describes the evolution of the [C/N] ra¬ 
tio during the EDU, provides a theoretical calibration of the 
[C/N]FDU-[M/H]-t relationship for low-mass stars which is com¬ 
pared with observations, and discusses several theoretical uncer¬ 
tainties relevant to the calibration. Section 3 compares our results 
with [C/N]FDU-[M/H]-t relationships determined from indepen¬ 
dent sets of stellar models, and this is followed by a summary of 
our results. 


* As is customary, we follow the spectroscopic notation for the ra¬ 
tio [A/B] between the number fractions of elements A and B, given as 
[A/B]=log(nA/nB).iMr “ log(uA/uB)o. The value of [A/B] is the same if 
element mass fractions (usually employed in stellar evolution calcula¬ 
tions) are used instead of number fractions (typically determined from 
spectroscopic observations). 
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Fig. 1. Upper panel: run of the [C/N] abundance ratio within 
a 0.9Mo stellar model ([M/H]=-0.35, scaled solar metal mix¬ 
ture), before the occurrence of the FDU (thick line) and right 
after the completion of the FDU (thin line). The horizontal axis 
displays the local fractional mass (mass enclosed within a dis¬ 
tance r from the centre divided by the total mass). Lower panel: 
as upper panel but for the individual mass fractions of N (dashed 
lines) and C (solid lines). 


2. Theoretical calibration of [C/N]fdu vs age 

Stellar model calculations show that at the end of the main se¬ 
quence (MS) phase, the outer convection zone progressively en¬ 
gulfs deeper regions, dredging to the surface matter which has 
been partially processed by H-burning during the MS (the FDU). 
Chemical elements affected by the FDU are essentially C, N, 
Li, and He. Both He and N abundances increase, whilst Li and 
C d ecrease, the *^C/^^C ratio d rops to values of ~20-25 (see, 
e.g.. lKarakas & Lattanzioll2Ql4 and references therein) and the 
i4n/15n ratio increases, due to mixing of the envelope’s initial 
composition with inner layers processed by the CN-cycle, where 
Li has also been burned. This slow variation in the photospheric 
chemical composition halts at the maximum penetration of the 
convective envelope, at which point the inner convective bound¬ 
ary starts to move back towards the surface, leaving behind a 
chemical discontinuity. This discontinuity, when crossed by the 
advancing (in mass) H-burning shell, causes the appearance of 
the RGB bum p in the luminosity func tion of old stellar popula¬ 
tions (see, e.g. lCassisi & Salarisll2013h . 

As an example. Fig. [T] displays the inner C and N stratifica¬ 
tions together with the corresponding [C/N] ratio before the start, 
and at the completion of the FDU respectively, of a Q.9Mq model 
with metal licity [M/H]--0.35 from the B aSTI stellar evolution 
databas43 (iPietrinferni et alJ |2004 l2006l) This value of [M/H] 
is obtained from Z=0.008, 7=0.256 (scaled solar metal mix¬ 
ture) according to [M/H]=log(Z/2£')jfar - log(Z/2£')o where Z, 


^ http://www.oa-teramo.inaf.it/BASTI/ 


X — 1 - Z - T, and Y denote the metal, hydrogen and He mass 
fractional 

The pre-FDU nitrogen chemical profile displays clearly three 
steps. Moving inwards, we find the initial abundance first, then 
the higher CN-cycle equilibrium value, and the even higher full 
CNO-cycle equilibrium value in the very central layers. The car¬ 
bon abundances follow their own pattern of initial, CN-cycle 
and CNO-cycle equilibrium abundances. After the FDU is com¬ 
pleted, convection has mixed the surface with the CN equilib¬ 
rium layers, and the final uniform abundance of N in the enve¬ 
lope has increased, whilst the C abundance has decreased. The 
corresponding behaviour of the [C/N] ratio is displayed in the 
upper panel of the same figure. 

At FDU completion the convective envelope contains ~72% 
of the total mass. When the mass of the model increases, surface 
convection at the FDU engulfs a larger fraction of the total stellar 
mass, and the net result is a decrease in [C/N] with increasing 
RGB mass -hence decreasing age of the host population- in the 
relevant mass (age) range. 

Figure|2]displays our theoretical calibration of the [C/N]fdu- 
[M/H]-t relation, obtained from the BaSTl stellar model library, 
including overshooting from MS co nvective cores and a sca led- 
solar heavy element distribution (see IPietrinferni et'ani2004 for 
more details). The MS convective cores have been extended be¬ 
yond the boundary given by the Schwarzschild criterion by an 
amount A„^.Hp, where Hp is the pressure scale height at the 
Schwarzschild border and A,„, a free parameter. For masses larger 
than or equal to 1.7Mq, BaSTI models employ Aov-Q.20Hp-, 
for stars less massive than LIMq /lo,,=0, while in the inter¬ 
mediate range of the model grid (M=l.l, 1.2, 1.3, 1.4, 1.5 
and 1.6Mn, respective ly) do,, yaries as Aov-{MjM q-Q.9)IA- (see 
IPietrinferni et aljr2004 for details). 

To determine the [C/N]FDU-[M/H]-t calibration we haye 
considered all low-mass model calculations at a giyen [M/H], 
and taken age and surface C and N abundances predicted for 
eac h mass at the RGB bump. The solar (C/N) ratio comes from 
the iGreyesse & Noelsl (1 19931) solar metal mixture employed in 
BaSTI calculations. Notice that the age range between FDU 
completion and RGB bump is negligible. 

For t aboye ~A Gyr there is a clear dependence of [C/N]fdu 
on [M/H], that increases with incresing age. Lower metal- 
licities imply higher [C/N]fdu^ due to the shallower conyec- 
tiye enyelopes and a consequent smaller oyerlap with the CN- 
equilibrium region. At any giyen metallicity the dependence of 
[C/N]fdu on t is highly non linear, the slope A[C/N]FDu/At in¬ 
creasing sharply for ages below ~5 Gyr; at any age this slope 
is shallower with increasing metallicity. As a consequence, the 
yalue of [C/N]fdu becomes increasingly less affected by [M/H] 
when t decreases, and it is predicted to be almost independent of 
[M/H] when t is below ~3 Gyr. 

[C/N]fdu measurements alone cannot howeyer proyide yery 
precise ages of indiyidual field stars, w hen considering typi¬ 
cal errors of ~0.10-0 .15 dex (see, e.g.. iBragaglia et ^ 2014 
iMikolaitis et^ 120111) . At [M/H]=-2.27 [C/N]fdu varies by 
~0.60 dex for a total age range equal to ~13 Gyr, whilst at 
[M/H]= 0.26 the [C/N]fdu yariation oyer the same age range is 
just ~0.35 dex. On the other hand, [C/N]fdu can proyide addi¬ 
tional age constraints on indiyidual RGB stars to be added to in- 


^ In general, if the metal mixture of the models is 
scaled solar, [M/H] = [Fe/H]. In case of a-element enhanced 
([a/Fe]>0) metal distributions, to a good approximation 
[M/H] ~ [Fe/H] -t log(0.638U -t 0.362) where L = 
dSalaris et all 199.31) . 
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formation coming from the log(g)-log(7’e//) diagram (and even¬ 
tually masses determined from asteroseismology), and also a 
general chronology for the formation of composite populations, 
when la rge samples of rC/NIppyi mea surements are available, as 
done bv iMasseron & Gilmorel (l2015h . We have also considered 
the [C/ NIpnii values predicted b y the BaSTI a-enhanced stellar 
models dPielrinferni et al.|[2006l) . and found a [C/N]FDu-[M/H]-t 
relation almost identical to the scaled solar one. 

The following analytical fit (our reference calibration) pro¬ 
vides the age of a star with measured [C/N]fdu [M/H], as 
obtained from the BaSTI models discussed above: 


t(Gyr) = 34.99 -h 156.3[C/N]fdu + 20.52[M/H]-h 
298.6([C/N]2pu) + 78.16 [C/N]fdu[M/H] + 7.82([M/H]2) 
+305.9([C/N]3nu) + 85.02([C/N]2pu)[M/H] 
-h22.93[C/N]fdu([M/H]2) -h 0.8987([M/H]^) 

+ 141.1([C/N]4pu) + 21.96([C/N]3pu)[M/H] 

+ 16.14([C/N]2pu)([M/H]2) + 1.447 [C/N]fdu([M/H]3) 

with an accuracy within +15%, for ages between 0.9 and 
14.5 Gyr, and the [M/H] range of Fig.|2] It can be applied to stars 
fainter than the RGB bump, which have completed the FDU. 
The values of surface gravity as a function of Teff at the comple¬ 
tion of the FDU (gpDu) and the RGB bump (gbump) predicted by 
the BaSTI models are given by: log(gFDu) = -11.951og(7’e//) + 
47.10 and log(gbump) = -14.441og(7’e//) +55.30, with ~0.15 dex 
spread, for the [M/H] range of Fig|2] Observed RGB stars that in 
a g-Teff diagram lie between these relationships for gpou and 
gbump, can therefore be age-dated with our [C/N]FDu-[M/H]-t 
calibration. 

The models employed for our reference [C/N]fdu-[M/H]- 
t relation do not include the effect of atomic diffusion and ra¬ 
diative levitation at old ages. Observations suggest that these 
processes are largely inhibited -at least from the convective 
envelopes- in low mass met al poor stars (see, e.g., iKorn et ^ 
l2006tlMucciarelli et al.l201ll and references therein). Moreover, 
after the FDU the surface C/N ratio (and to a good approxima¬ 
tion all chemical abundances) is basically the same as for the 
case of neglecting dif fusion in low mass me tal poor models, as 
shown for example bv iMichaud et al.l (l2007h . The effect of fully 
efficient diffusion would be just a decrease by about 5% of the 
estimated ages at fixed [C/N]fdu and [M/H]. 

Our reference calibration does not account for overshooting 
from the lower boundary of convective envelopes. To check the 
impact of this process on [C/N]fduj we have computed selected 
additional models at [M/H] =0.06 (Z=0.02, 7=0.273) extending 
the fully mixed envelope by a distance To,, = Q.3Hp, a choice 
consistent with indications from the discrepancy between pre¬ 
dicted and observed brig htness of the RGB bump in Galactic 
globular clusters (see, e.g. lCassisi et alJ 201 ih . We found that the 
change (decrease) of [C/N]fdu is generally very small, equal at 
most to ~ 0.02 dex. 

We performed three additional tests at [M/H] =0.06. The first 
one is an assessment of the effect of neglecting convective core 
overshooting during the MS. Small differences in [C/N]fdu at a 
given t appear for t below ~3 Gyr, and reach at most 0.05 dex 
([C/N]FDU higher in models without overshooting) at t=l Gyr. 
We have then calculated models to study the effect of varying 7 
(a decrease AY = 0.02) and the mixing length parameter (a 
decrease in a^x by 0.3 compared to the solar calibrated reference 



Fig. 2 . Theoretical predictions for the run of [C/N]fdu as a func¬ 
tion of age for selected metallicities. From top to bottom, at an 
age of 10 Gyr, [M/H] is equal to -2.27 (dotted line), -1.49 (short 
dashed), -1.27(long dashed), -0.66 (dot-dashed), -0.35 (thick 
solid), 0.06 (thin solid) and 0.26 dex (thin long dashed), respec¬ 
tively. Data for halo field stars and the labelled open clusters are 
also plotted (see text for details). 


value ami_o=2.01 of the BaSTI calculations) respectively. In both 
cases the predicted [C/N]FDU-[M/H]-t relation is unchanged. 

We have then compared our reference theoretical calibration 
with observations of [C/N]fdu in field halo stars and a few open 
clusters of different ages. Observational data are scarce and error 
bars sizable when compared to the scaling of [C/N]fdu with age 
and metallicity, and do not allow a stringent test of the models. 

To sample old and meta l-poor system s we c onsidered the 
field halo stars observed by iGratton et al.l (l200r)h . with an es¬ 
timated mean [C/N]fdu = -0.07 + 0.05 (error equal to the Icr 
dispersion around the mean), assuming an age of 13.5 Gyr for all 
the stars. These objects span the [Fe/H] range between ~ -1.3 

and-2.0 dex -corresponding to [M/H] b etween ~ -1.0 and ~ 

-1.7 dex for a typical [a/Fe]~0.4, (see e.g. lCarretta et alj|2000l 
and references therein)- a regime where theory predicts a very 
small variation in [C/N]fdu with metallicity, even at old ages 
(see Fig. | 2 ]i hence it is justified to consider a mean [C/N]fdu 
for the sample. We t hen co nsidered three stars in NGC6791 
(iBragagliaet a .112014 . M67 dGilrov & Brownlll991h . NGC6819 


(ICarlberg et al. 


12015l excluding the Li-rich star that appears on 


the lo wer RGB of the cluster), and NGC752 ([B 6 cekTo 2 cue^^ 


201^) respectively, plus one star in NGC2506 (iMikolaitis et al 


201 11) . In case of multiple objects from the same cluster we took 


the mean value of [C/N]fdu and the dispersion around the mean. 
The evolutionary status of the observed open cluster stars (post 
FDU but pre-RGB bump) has been determined making use of 
the appropriate isochrones from iPietrinferni et al.l (1200411 . 

F or NGC6791 w e empl oyed the distance modulus and age 
from iBrogaard et aP (l2012h . and [Fe/H] =0.37 dBragaglia et al.l 
l2014t) . In case of M67 we considere d [Fe/H] =0. 05, dis¬ 
tance modulus and age all taken from iBellini et al.l (l2010l) . 
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whilst for NGC6819 we adopted [Fe/H]=0.09, distance and 
age from ISala ris et al. (|2009|). and for NGC752 we employed 
rFe/H1^-0.02 (lB6cekToDcuetal.ll2015h, t^l.7+0 Gyr, ( w - 


M)v=8.36 (see discussion in iBocek Toncu et al.l [20151). In 


=-0.24 d Mikolait is et al.l 

Marconi et al.l (1199^ 


case o f NGC2506 we assumed [Fe/H 
UMl, age and distance modulus from 

For M67, NGC6819 and NGC2506, where the original pa¬ 
pers either determined only the number ratio (C/N) (the case 
of M67 and NGC6819) or assumed (not determined with the 
same methods) a value for (C/N )0 from independent estimates 
of the solar (C/N) rat io (the case of NGC2506 ). we employed 
log(C/N)o=0.58 from iGrevesse & Noelsl(ll993l) . as in the BaSTl 
calculations. This causes an increase in the measured [C/N] 
by just 0.02 dex for NGC250 6. co mpared to the value pro¬ 
vided by Mikolait is et alJ (1201 ll) . Both lAsplund et akl (l2009l) and 
iGrevesse & Sauva |(ll998l) solar mixtures give log(C/N)o=0.60, 
whilst icaffau et all(l201 111 estimate gives log(C/N)o=0.64. 

Figure |2] summarizes the comparison between the observa¬ 
tional data described above and our reference theoretical predic¬ 
tions. As anticipated before, the data cannot provide a stringent 
test of the theory, but overall they confirm the predicted decrease 
in [C/N]fdu with age. 


3. Comparison with other steiiar modeis 

We compared our referen ce [C/Nlpnii-[M/Hl -t calibration with 
results from the PAR SEC dBressan et alJ201^ and STARE VOL 
(iLagarde et alj|201^ databases of stellar evolution calculations, 
which provide also surface chemical abundances of several el¬ 
ements during the evolution of models of different mass and 
initial chemical composition. The solar heavy element distribu¬ 
tion differs amongst BASTI, PARSEC and STAREVOL. BaSTI 
calculations employ the IGrevesse & Noelsl dl993h distribution, 
whilst PARSEC and STAREVOL models adopt the lCaffau et alJ 
d2Ql ill and lAsplund et al.l d2009t) determinations, respectively. 
These three different estimates provide different values of the 
solar metal mass fraction, that we accounted for in the calcula¬ 
tions of [M/H], whilst the solar C/N ratio -which is the C/N ratio 
of the models’ initial chemical composition and is also needed 
to calculate [C/N]fdu- does not change by more than 0.06 dex, 
as discussed in the previous section. PARSEC and STAREVOL 
calculations also include convective core overshooting during 
the MS as our reference calibration, although with different ef¬ 
ficiencies. In addition, PARSEC models account for overshoot¬ 
ing from the convective envelope and atomic diffusion, partially 
inhibited from the surface convective layers, as discussed in 
iBressan et al.l (l2012h . 

The top panel of Eig. [3] compares our [C/Nlppu-t relations 
for [M/H] =0.06 (Z=0.02, Y=0.273) and -0.66 dex (Z=0.004, 
Y=0.251), with PARSEC results for [M/H]=-0.02 (Z=0.014, 
Y=0.273) and -0.58 dex (Z=0.004, y=0.256) - the BaSTI and 
PARSEC models are not calculated for exactly the same [M/H] 
grid. There is good quantitative agreement between the two sets 
of results for ages larger than ~ 9 - 10 Gyr, that deteriorates to¬ 
wards younger ages. In general PARSEC models predict lower 
[C/N]fdu values at a given t and [M/H], typically by ~0.05 dex 
at t~5 Gyr, increasing to ~0.10 dex at t~1.0 Gyr. However, 
what is more relevant is the corresponding age difference for 
a given [C/N]fdu^ at fixed [M/H]. PARSEC ages are systemati¬ 
cally higher, with differences of the order of ~ 1.5 - 2.0 Gyr for 
t between ~4 and ~8 Gyr. At older ages the difference decreases 
sharply to almost zero, whereas it is of the order of ~0.5-l Gyr 
at ages below ~4 Gyr. 


By comparing [C/N]fdu and EDU ages for various masses at 
fixed [M/H], we found that the discrepancy is mainly driven by 
lower [C/N]fdu values at fixed mass for the PARSEC models. 
As a test we then calculated some models with MS core over¬ 
shooting (without diffusion nor envelope overshooting) as in our 
BaSTI reference calibration, but with the lCaffau et alJ(l201 Ih so¬ 
lar metal distifbution, for Z=0.0172 and T=0.270, which is the 
initial solar composition determined by calibrating the standard 
solar model (with diffusion included). This initial chemical com¬ 
position corresponds to [M/H] =0.04, and the resulting [C/N]fdu 
values as a function of age are also reported in Eig. [3] Just the 
effect of the solar metal distribution does not seem to be able 
to explain the differences with PARSEC results at the similar 
metallicity [M/H]=-0.02. 

PARSEC calculations include overshooting from the con¬ 
vective envelopes that increases when the s tellar mass increases 
(see the discussion in IBressan et alJl 2012 l) . and we tested also 
this effect. By employing the same envelope overshooting pre¬ 
scription as in the PARSEC calculations -Tov = 0.7 x Hp for 
masses above 1.4Mo, decreasing at lower masses down to 
Tov = 0.05 X Hp when M< I.IMq)- we obtained a decrease in 
[C/N]fdu at fixed age of the order of 0.03 dex for ages below 
~3Gyr, this variation getting smaller at higher ages. Adding the 
0.03 dex decrease in [C/N]fdu to the square symbols in Eig.[3]at 
the ages ~1 and ~3 Gyr moves the points closer to the PARSEC 
results, still without matching them. In conclusion, part of the 
disagreement between the two calibrations seem to be due to the 
different metal distribution and the effect of envelope overshoot¬ 
ing, even though some additional structural differences in the 
models -maybe related to the different efficiency of convective 
core overshooting during the MS- play a role. 

The lower panel of Eig. [3 displays a comparison with the 
STAREVOL calculations for [M/H]=0.0 (Z=0.014, Y=0.266) 
and [M/H]=-0.56 (Z=0.004, Y=0.253), respectively, without 
and with the inclusion of rotation. Both sets of calculations 
include also thermoaline mixing, that sets in after the RGB 
bump, and is irrelevant for our comparisons. The shape of the 
[C/N]FDU-t relation is less smooth than the case of BaSTI mod¬ 
els, due to a coarser grid of masses. STAREVOL calculations 
without rotation predict lower [C/N]fdu at fixed age and [M/H], 
typically by ~0.05 dex. The main culprit for these differences is 
the longer lifetime of STAREVOL models of a given mass and 
[M/H]. Eor a given [C/N]fdu, at fixed [M/H]. STAREVOL cal¬ 
ibration provides therefore ages systematically higher than our 
calibration, with differences of the order of ~1.5 Gyr, decreas¬ 
ing to ~ 0.5 for ages below ~2 Gyr. 

The STAREVOL models with rotation have been calculated 
with an initial rotational velocity equal to 45% of the critical 
velocity for each mass, and display much lower [C/N]fdu val¬ 
ues at a given t and [M/H], with differences up to ~0.20 dex 
compared to our reference calibration. This is due mainly to ro¬ 
tational mixing that pushes N from the inner regions at CNO 
equilibrium (higher N abundance at CNO equilibrium compared 
to the NO equilibrium) into overlying layers later to be engulfed 
by convection during the EDU. At the same time C is pushed 
inwards from the envelope, given that CN and CNO equilibrium 
abundances are lowe r than the initial ones (see, e.g., Fig.l 1 from 
iLagarde et al.ll201^ . The effect of rotation on lifetimes at the 
EDU is very small. In this mass range rotation causes a small in¬ 
crease in the MS lifetime, tog ether with a shorteni ng of the RGB 
timescales (see discussion in ILagarde et al.ll 2012 h . compared to 
non-rotating models. 

A comparison of Eigs. [3]and shows that at the lower age 
end of the [C/N]FDU-t calibration the existing sparse observa- 
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Fig. 3. Comparison between the BaSTI (thick lines) and 
PARSEC (thin lines, upper panel), STAREVOL (thin lines 
lower panel) predictions for the [C/Nlpou-t relation at selected 
metallicities. BaSTl results for [M/H]=0.06 (solid line) and 
[M/H]=-0.66 (dash-dotted line) are displayed together with 
PARSEC results for [M/H]=-0.02 (solid line) and [M/H]=-0.58 
(dash-dotted line), and STAREVOL results for [M/H]=0.0 (solid 
line) and [M/H]=-0.56 (dash-dotted line). STAREVOL models 
with rotation for the same metallicities are also displayed (with 
the same line styles). In the upper pa nel squares denote B aSTI 
calculations for [M/H]=0.04 and the ICaffau et alJ (l201 ih solar 
heavy element distribution (see text for details). 

tional data cannot strongly confirm or disprove any of these al¬ 
ternative calibrations. However, at old ages models with rotation 
seem to largely underestimate [C/N]fdu compared to NGC6791, 
if we assume that the scaling of [C/N]fdu with [M/H] is simi¬ 
lar to BaSTI results between solar and twice solar metallicity 
(almost negligible difference), given that [M/H]=0 is the upper 
boundary of STAREVOL metallicity range. 

4. Summary 

We have analysed the use of [C/N]fdu as age indicator for RGB 
stars, and derived a theoretical calibration of [C/N]fdu in terms 
of total metallicity and age, for ages larger than t ~1 Gyr, based 
on the BaSTI stellar models with convective core overshooting. 
In terms of general qualitative properties, lower metallicities im¬ 
ply higher [C/N]fdu at old ages, but [C/N]fdu becomes increas¬ 
ingly less affected by [M/H] when t decreases, being almost in¬ 
dependent of metallicity below ~3 Gyr. 

We have tested the robustness of our calibration against vari¬ 
ations in the initial heavy element distribution (scaled-solar vs a- 
enhanced), efficiency of overshooting from MS convective cores 
and from the convective envelopes, variations in the initial He 
abundance and in the mixing length parameter am\- The small 
sample of available measurements of [C/N]fdu at least qualita¬ 
tively confirm the overall trend of the predicted [C/N]fdu with 
age and [M/H]. 


Eor ages ranging between 1 and 14 Gyr, [C/N]fdu changes 
by at most ~0.60 dex at the lowest metallicities, and by just 
~0.35 dex at supersolar [M/H]. Eor typical uncertainties of 0.10- 
0.15 dex on the observed abundance ratios, [C/N]FDU-[M/H]-t 
relationships cannot therefore give very precise ages of individ¬ 
ual stars, but may provide a general chronology for the forma¬ 
tion of composite populations, when large samples of RGB stars 
with [C/N]fdu measurements are available, and add further con¬ 
straints to analyses of RGB samples from the log(g)-\og(Teff) 
diagrams. 

Comparisons with the [C/N]FDU-[M/H]-t relations from 
the PARSEC stellar models display [C/N]fdu differences at 
fixed t and [M/H] that increase with decreasing age, whilst 
STAREVOL models without rotation show a more constant 
offset of [C/N]fdu- Age differences up to ~2 Gyr (PARSEC 
and STAREVOL ages being systematically higher) for a given 
[C/N]fdu [M/H] pair are possible, mainly in the age range 
between ~3-4 and ~8 Gyr. STAREVOL calculations with ro¬ 
tation display a further decrease in [C/N]fdu at a given t and 
[M/H], that has an even larger impact on the ages derived from 
[C/N]fdu. 

Additional accurate spectroscopic measurements of 
[C/N]fdu in star clusters with well established ages and metal¬ 
licities are necessary to better test the theoretical calibrations of 
this age indicator. 
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